Intruders in the Dust: Air-Driven Granular Size Separation 
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Using MRI and high-speed video we investigate the motion of a large intruder particle inside a 
vertically shaken bed of smaller particles. We find a pronounced, non-monotonic density dependence, 
with both light and heavy intruders moving faster than those whose density is approximately that 
of the granular bed. For light intruders, we furthermore observe either rising or sinking behavior, 
depending on intruder starting height, boundary condition and interstitial gas pressure. We map 
out the phase boundary delineating the rising and sinking regimes. A simple model can account for 
much of the observed behavior and show how the two regimes are connected by considering pressure 
gradients across the granular bed during a shaking cycle. 

PACS numbers: 45.70.Mg, 64.75.+g, 83.80.Fg 



Unlike thermal systems which favor mixing to increase 
entropy, granular systems tend to separate under an ex- 
ternal driving mechanism such as vibrations Q, Q . This 
is commonly known as the Brazil Nut Effect, in which 
a large particle, the "intruder", rises to the top of a 
bed of smaller background particles 0, 0, More re- 
cently, new behavior was discovered for the limit of very 
small bed particles ( "dust" ) , in particular the sinking of 
light intruders [|| , and a non-monotonic dependence 
of the rise time on density 0, 0, El • A number of the- 
ory and experimental papers ex plor ed different aspects of 
this surprising behavior [ll|, |l2|, yjj LUj , but so far- 
there has been no consensus about either the underlying 
mechanisms or the relative importance of various system 
parameters in driving the intruder motion. 

Here we present results from a systematic investigation 
of both the intruder motion and the bed particle flow. 
Our central finding is that there is a phase diagram which 
delineates rising and sinking behavior of the intruder as 
a function of interstitial gas pressure, intruder density, 
and initial intruder height within the container. Our re- 
sults lead to a physical model that provides a unifying 
framework to describe both rising and sinking regimes. 
In this way, the work presented here connects previously 
disjointed pieces of a puzzle that pointed to the impor- 
tance of pressure gradients |3, la li| but approached the 
two regimes as separate phenomena. As a consequence, 
our findings directly contrast with the mechanisms pro- 
posed in Refs. |g, |lfj, [lj, |l6( that neglect interstitial gas 
flow. 

We placed granular material inside an acrylic cylinder 
(inner diameter 8.2 cm) mounted on a shaker and used 
individual, well-spaced sine wave cycles ("taps") of fre- 
quency / and amplitude A to vibrate the vessel vertically. 
The cell could be evacuated to a gas pressure, P. Both 
smooth and rough cells (created by gluing glass beads 
to the interior walls of an otherwise smooth cell) were 
used to study the effect of wall friction. A large intruder 



sphere of diameter D was buried in the bed of background 
spheres (diameter d) at a height h s measured from the 
vessel bottom to the intruder top (See Fig.^a)). A range 
of diameter ratios D/d, shaking parameters, and back- 
ground particle materials (glass, zirconium oxide, tapi- 
oca, and seeds [for MRI]) were investigated. The intruder 
density, p, could be tuned. We measured the number of 
taps required for the intruder to break through the up- 
per free surface, T r i SC , or to reach the bottom, T s j n k. In 
addition, the intruder position could be recorded with 
high-speed video throughout the shaking process by at- 
taching to the intruder a thin, vertical straw extending 
above the upper surface. We verified that the straw did 
not affect the intruder motion. 




FIG. 1: MRI image of a large intruder sphere rising in a bed 
of smaller particles. Images are vertical cut through the center 
of the cylindric cell, (a) Layered bed and the starting position 
of the intruder before shaking, (b)-(d) Size separation in the 
presence of convection for intruders of three different densities 
at r = 5, / = 13 Hz: (b) p/p m = 0.08 after 23 taps; (c) 
P/Pm = 0-44 ~ P*/p m after 30 taps; and (d) p/p m = 2.33 
after 11 taps. Here p m = 1.10 g/ml is the average density of 
the seeds used as bed particles. 



To track the evolution of both the intruder and the sur- 
rounding bed material we used magnetic resonance imag- 
ing (MRI). Fig.njshows for the rough cell how the motion 
of the intruder is coupled to that of the background con- 
vection. Horizontal layers of poppy seeds (MRI active) 
were alternated with layers of rajagara seeds (MRI in- 
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active), with the intruder initially placed on the lowest 
black layer (Fig. ^a)). Figures ^b)-(d) show the situ- 
ation, for intruders with progressively higher densities, 
after several taps when the intruder had risen slightly 
more than its own diameter. In both Fig. ^b) and^d), 
the intruders move faster relative to the dark layer on 
which they were originally placed. In Fig.^c), the inter- 
mediate weight intruder with density close to the effect 
density of the bed remains at the same relative position 
to that layer. Thus, its rise time corresponds to the con- 
vection speed of the bed particles immediately below it. 
Above the intruder in Figs, ^b) and^c), the alternat- 
ing light and dark layers are no longer horizontal, because 
convection speeds up near the free surface and con- 
vection rolls churn the material into a swirl pattern near 
the side walls. In Fig. ^d), the dense intruder rises so 
rapidly that these convection rolls have not had a chance 
to move the background material appreciably. As the 
dense intruder rises, it pushes a wedge-shaped volume of 
material above it and creates a wake below it. In Fig.^b) 
for the light intruder, there is no sign of a similar wake. 

One aspect common to prior observations of the sink- 
ing regime is the need for small bed particle sizes, typi- 
cally d < 0.5 mm However, this is precisely also 
the size regime in which it becomes very difficult to ex- 
cite axi-symmetric bed flow patterns yjj- In general, 
wall-driven convection sets up axi-symmetric convection 
rolls in the bed However, when either bed perme- 
ability or wall-driven friction is sufficiently reduced, the 
convection becomes asymmetric 0, so that intrud- 
ers are driven towards the cell wall. We find the onset 
of asymmetric flow for d < 0.35 mm in the rough cell 
and d < 0.5 mm in the smooth cell. In the following, 
we focus on representative results obtained with d = 0.5 
mm glass beads (p m = 2.5 g/ml), D = 25 mm intruders, 
bed fill height 85 mm, and dimensionless shaking accel- 
eration T = A(2irf) 2 /g = 5 at / = 13 Hz. The two cell 
types then allow us to investigate how the overall bed 
flow affects the intruder motion. 

Despite the clearly different flow patterns in the two 
cells, we find a qualitatively similar dependence of the 
intruder rise time, T r i se , on the ratio of intruder den- 
sity, p, to the bed material density, p m (Fig. EJ. For 
both boundary conditions, at ambient pressure, T r ise has 
a pronounced peak at p* ~ 0.5p m . On either side of 
this peak, there is a large increase of the intruder veloc- 
ity (i.e., a decrease in T riso ). When P is decreased, the 
peak amplitude shrinks and p* moves to smaller values 
as shown for both cells in the inset to Fig. |3[a). From 
MM measurements (Fig. QJ, we know that at the peak 
p* the intruder rises with the convection, while intruders 
on either side of the peak rise faster. This result is in 
contrast with previous experiments 0, ^J, where con- 
vection is measured without the intruder. This indicates 
the presence of the intruder slows down the convection. 
When P is lowered, the convection speeds up so that 
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FIG. 2: (a) Intruder rise time T r i sc versus p/p m at differ- 
ent pressures P. (a) Rough cell: (■) 101 kPa, (•) 47 kPa, 
(▼) 13 kPa, (a) 6.7 kPa, (0) 0.13 kPa. Inset: rise time 
peak/divergence position p* / p m as a function of P in the 
rough (o) and smooth (■) cell, (b) Smooth cell: (■) 101 
kPa, (o) 27 kPa, (♦) 0.13 kPa. Inset: The sinking regime at 
2.7 kPa. In both cells, h s — 5.5 cm. 



the amplitude of the rise-time peak shrinks as shown in 
Fig. |21 For the rough cell, p* decreases more strongly 
than in the smooth cell and at low P it is indistinguish- 
able from p* = 0, having reached our lowest measurable 
density. At sufficiently low pressure (P < 0.13 kPa), the 
non-monotonic behavior in T r ; sc disappears and the curve 
is featureless. The slow increase in T r j se as p increases we 
believe is due to the heavy intruder burrowing back into 
the bed at the end of each cycle. Aside from this there 
is essentially no density dependence to the rise time. For 
p < p* and P ~ 2.7 kPa a dramatic change in behavior 
can be observed, seen most clearly in the smooth cell. As 
shown in the inset to Fig.|2Jb), at P = 2.7 kPa, instead of 
speeding up again as p is lowered below p*, the intruder 
stops rising and begins to sink. In this case, the peak 
in T r ; sc at p* turns into a discontinuity as shown in the 
inset of Fig.^b). We find similar behavior in the rough 
cell but over a much smaller region of parameter space. 
We find that the transition between the two regimes of 



3 




FIG. 3: Phase diagrams delineating the rising and sinking 
regimes for d = 0.5 mm beds at various pressures P. Above 
each phase boundary, intruders rise. Shaded area shows the 
sinking regime at ambient pressure. Main panel: (■) 101 
kPa, (x) 27 kPa, (a) 2.7 kPa, (V) 0.67 kPa in the smooth 
cell. Inset: (■) 101 kPa, (o) 27 kPa in the rough cell. In 
both cells, all intruders rise from the bottom (h s — D) when 
P < 0.13 kPa. The dashed lines indicate h s of Fig.H 



intruder rising and sinking is controlled by at least five 
main parameters: p/p m , P, d, h s and cell wall roughness. 
This allows us to construct phase diagrams for different 
wall roughness as shown in Fig. where we plot the 
maximum starting height for sinking, h c , as function of 
relative density. Results for different d will be presented 
elsewhere If the intruder is placed close enough 

to the top surface, it will invariably rise. The regions 
where the intruder sinks depend on pressure, P, and, 
when they occur, happen at low relative densities, pj p m , 
and low heights, h s . Interestingly, the sinking regime is 
largest for intermediate pressures around 1 — 30 kPa. At 
lower pressure, this regime disappears completely and all 
particles rise at all depths as they are swept along via the 
background convection rolls. We find the sinking regime 
only for p < p* . 

In order to determine how the intruder interacts with 
the background, we measured its motion during flight 
using high-speed video. In Fig. 0] we plot the vertical 
displacement Ah (measured in the moving frame of the 
cell) versus time. The inset of Fig. 0Ja) shows intruder 
trajectories in the lab frame. In both rough and smooth 
cell, the bed lifts off the cell bottom at t = 0, the particles 
near the walls land at 0.05 s and the condensation front of 
the bed reaches the intruder at ~ 0.07 s. Between liftoff 
and landing of the bed, the heavy intruder in the rough 
cell follows a parabolic trajectory with a total downward 
acceleration only slightly greater than g (Fig. 0] inset) 
whereas, as is evident from the larger effective downward 
acceleration, the lighter intruder experiences larger drag 
due to the air flow. Before landing, the light particles 
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FIG. 4: High-speed video intruder displacement in the mov- 
ing frame. The vertical dashed lines delineate the part I and 
part II of the period (see text), (a) Three different densi- 
ties with h s = 7.0 cm in the rough cell. (■) p/p m = 0.043, 
(□) p/p m = 0.52 ~ p*/p m , (A) p/pm = 3.3. Inset: Same 
trajectories in the lab frame. The "kink" for the lightest in- 
truder is highlighted by the dashed circle, (b) Light intruder 
[p/pm = 0.043) in smooth cell rising at h 3 = 6.5 cm > h c (■) 
and sinking at h s — 3.5 cm < h c (•). The third curve (a) 
tracks the gap underneath the bed. 



show a sharp change in behavior. This is seen as a kink 
in the curves occurring near 0.045 s (apparent in the lab 
frame). This kink, which slows down the landing of the 
light particles, appears to be responsible for the increase 
in upward velocity for p < p*. In the smooth cell, the 
kink is not as dramatic as in the rough cell but the de- 
parture from the downward parabolic profile is still ev- 
ident and the intruder clearly moves down slower than 
the bed. When a light intruder is in the sinking region 
of the phase diagram, its trajectory, as seen in Fig.^Jb), 
is pulled down with respect to the bed in the first part 
of the shaking period. 

These results clearly identify the interstitial gas as the 
cause of the non-monotonic behavior in T r i SC and of the 
reversal from rising to sinking. They contradict several 
recent models that ignore air effects or treat the ri sing 
and sinking regimes as unrelated phenomena [lOj, Il41 |l6( . 
However, we can explain the rise/sink cross-over in the 
phase diagram and predict quantitatively the peak at 
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T r ise{p* /Pm) with a model that treats the bed as a porous 
piston with permeability k. We consider two parts of 
the free-flight portion during each shaking cycle (vertical 
lines in Fig. 0J: In part I, gas flows down into the gap 
opened up at the cell bottom and produces a drag that 
adds to the gravitational acceleration g. In part II, the 
gap closes, so that the gas pushed upward through the 
bed creates a drag force opposing the inertial force. 

We approximate both the packing fraction / of bed 
particles and the pressure gradient to be constant 
throughout. A more detailed discussion will be presented 
elsewhere [19J. By Darcy's law, incompressible gas of 
viscosity p and velocity u flowing through the bed pro- 
duces a pressure gradient dP/dz — pu/k. This gradient 
leads to a drag force F4 = (pu/k)V on a bed volume V. 
Similarly, an intruder of volume Vi experiences a force 
§ l P{z)dS = {dP/dz)V = {fiu/k)Vi. Consider the free 
flight motion of an intruder of mass = pV. relative to 
a neighboring, identical volume V. comprised only of bed 
material. Depending on whether the sign of the mass 
difference Am = (p — / p m )Vi between the two volumes is 
positive (negative), the intruder during part I will push 
against the material above (below) it. We calculate the 
net force on the compound object consisting of the in- 
truder and the vertical column (diameter D) of bed ma- 
terial either above or below it. We find that, during part 
I, intruders with p/p m < / will sink, while those with 
PI 1 pm > / will rise relative to the neighboring bed. 

During part II of the free flight the situation changes. 
Now the heavy (light) intruder accelerates together with 
the column of material below (above) it. Because in- 
ertial and drag forces oppose each other, light intruders 
experience smaller acceleration magnitudes and fall more 
slowly than the surrounding bed. In particular, near the 
top surface light intruders can reach their terminal veloc- 
ity before the bed collides with the base. We see evidence 
of this in the inset to Fig. Eta) as the parabolic trajec- 
tory changes to a linear, constant velocity segment in the 
lab frame. Finally, we assume that any gap around the 
intruder is immediately filled by bed particles so that its 
displacement over one shaking cycle is the sum of the 
displacements from parts I and II. 

From these considerations several predictions emerge. 
First, at p* / ' p m — f the intruder experiences no motion 
relative to the bed. This is what is seen in Fig. ^c) 
and in the smooth cell leads to a divergence in T r ; sc . In 
the presence of convection, this divergence should be cut 
off by the convective rise time. Second, the peak in T r i se 
should occur at a packing fraction value corresponding to 
a loosely packed "in-flight" bed configuration. Thus, we 
expect similar values for p* j p around / ~ 0.5 using dif- 
ferent bed particle sizes and shaking parameters |ill9lllo|. 
Third, for light intruders the amount of sinking during 
part I and rising during part II should depend on their 
vertical position in the bed. Thus, for pj p m < f there 
should be a critical initial height, /i c , separating rising 



and sinking behavior. Because convection will produce 
a bias towards rising, it will reduce h c in the rough cell. 
Qualitatively, this explains the key features of the phase 
diagrams in Fig. Finally, the model implies that the 
sinking found by us and others 0, depends on a pres- 
sure gradient and should vanish in vacuum, where only 
the convection-driven, density-independent rising effect 
survives This is confirmed by the data for P < 0.13 
kPa in Fig. El 

While the above model provides a mechanism for the 
demise of the density dependence with decreasing pres- 
sure it neglects a number of aspects that might be impor- 
tant: variation in dP/dz, compressibility of the bed and 
associated variations of / with position, and pressure- 
dependence of convection. Thus, the model cannot pre- 
dict the detailed shape of the phase diagram (Fig. [3| nor 
the inset of Fig. El^a). From earlier results |lj| by the 
Duke group, we would have expected a significant pres- 
sure dependence only once P approached values close to 
2 kPa. However, our results clearly demonstrate changes 
in the behavior at much larger P. 

By presenting the phase diagram as a function of 
various system parameters, we showed that the pre- 
vious perplexing results obtained by different groups 
H> IE El belong to different regimes of the same 

phenomenon. Using MRI and high-speed video, we found 
the underlying mechanisms of the density-dependent be- 
havior of intruders. Moreover, a simple model can give a 
qualitative description of the key experimental results. 
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